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Between May 2009 and May 2010, the IceCube neutrino detector recorded 32 billion of at- 
mospheric muons generated in air showers produced by cosmic rays in the TeV energy range. 
With such high statistics sample it is possible to observe, for the first time in the southern 
hemisphere, an energy dependence in the Galactic cosmic ray anisotropy up to a few hun- 
dred TeV. This study shows that the same large-scale anisotropy observed at median energies 
around 20 TeV is not present at 400 TeV; the anisotropy observed at 400 TeV shows substan- 
tial differences with respect to that at lower energy. In addition to the large-scale features 
observed at 20 TeV in the form of strong dipole and quadrupole moments, the data include 
several localized regions of excess and deficit on scales between 10° to 30° . The features ob- 
served at both large and small scale are statistically significant, but their origin is currently 
unknown. 



1 Introduction 

The IceCube detector, deployed around 2000 m below the surface of the South Polar ice sheet, 
is designed to detect upward-going neutrinos from astrophysical sources. However, it is also 
sensitive to downward-going muons produced in cosmic ray air showers. To penetrate the ice 
and trigger the detector, the muons must possess an energy of at least several hundred GeV, 
which means they are produced by primary cosmic rays with energies in excess of several TeV. 
At these energies it is expected that interactions of cosmic rays with Galactic magnetic fields 
should completely randomize their arrival directions. Nevertheless, at TeV energies, multiple 
observations of anisotropy in the arrival direction distribution of cosmic rays havejaeen obse rved 
on large and small angular scales by the Tibet AS7 array^, Super-Kamiokande ^ , Milagro^^Hl^ 
ARGO-YBJ^. These measurements were performed in the northern hemisphere. The first ob- 
servation of the large-scale cosmic ray anisotropy in the southern sky has been reported by the 
IceCube detector 6 in 2010. 

The observations at several TeV are quite similar across all experiments and between the north- 
ern and southern hemispheres. However, at high energies observations from several experiments 
have been really controversial. For example, the Tibet AS7 collaboration has published a null 
result at 300 Tevffl, while at 370 TeV the EAS-TOP collaboration has claimed evidence of an 
anisotropy^. 

In this proceedings we describe the energy-dependence of the large-scale anisotropy using 
IceCube data. We also describe a search of the southern sky for anisotropy on small and medium 
angular scales. 



2 The IceCube Detector and Data 



IceCube, completed in December 2010, is a km'^-size neutrino detector frozen into the glacial 
ice sheet at the geographic South Pole. High-energy particles passing through the detector emit 
Cherenkov radiation and their tracks are recorded by an array of 5160 Digital Optical Modules 
(DOMs)EI embedded in the ice. The DOMs are attached to 86 vertical cables, or strings, deployed 
at depths between 1450 m and 2450 m below the surface of the ice sheet. Between 2009 and 2010, 
59 strings had already been deployed and the detector was running in this smaller configuration 
(IC59). 

The main trigger used for physics analysis in IceCube is a simple majority trigger which requires 
coincidence of 8 or more DOMs hit in the deep ice within a 5 /is window. The trigger rate of 
downgoing muons in IC59 was 1.7 kHz, a factor of 10^ larger than the neutrino rate. Due to the 
high rate of muon events, the muon tracks are reconstructed online using a maximum likelihood 
reconstruction, and the results are stored in a Data Storage and Transfer (DST) format to satisfy 
bandwidth constraints at the South Pole. 

The analysis presented in this paper uses the DST data collected during IC59 physics runs 
between May 2009 and May 2010. The data set contains approximately 3.4 x 10^*^ muon events 
detected with an integrated live time of 334.5 days. A cut in zenith angle {9 > 115°) is used 
to remove misreconstructed tracks near the horizon, reducing the final data set to 3.2 x 10^*^ 
downgoing events. Simulations indicate that the median angular resolution is 3° and the median 
primary energy of the cosmic ray data set is 20 TeV. The primary cosmic ray energy resolution 
is of order 0.5 in Alog(£'). 



3 Data Analysis and Results 

3.1 Evolution of the Anisotropy with Energy 

Because we must measure cosmic rays indirectly via observations of muons produced in the 
extensive air showers, the energy of the cosmic ray primary particles is inferred using estimates 
of the muon energy. 

In this analysis the estimate of cosmic ray energy is based on the number of DOMs hit by 
Cherenkov photons (i.e. number of channels, or Nch)- Due to the zenith angle (6) dependence 
of the relation between Net and the cosmic ray primary energy, a two dimensional cut in Nch 
and 9 is used to divide the data into two energy samples. The low energy sample contains events 
with a median energy of 20 TeV, where 68% of the events are between 4 and 63 TeV; and the 
high energy sample contains events with a median energy of 400 TeV, with 68% of the events 
between 100 and 1258 TeV. 

We apply the same method used in the past analysis!^ to investigate the arrival direction dis- 
tribution of the cosmic rays at these two energies. Equatorial maps of relative intensity are 
produced for the two energy samples and shown in Fig. [TJ The sky is binned into an equal-area 
grid with a resolution of 0.9° using the publicly-available HEALPixf^^ library, and smoothed with 
3° to match the data angular resolution of the online muon reconstruction. 
To study the evolution of the anisotropy with energy, we project the maps in right ascension 
and fit the resulting ID distribution to a first and second-order harmonic function of the form: 
J2'j=i = Aj cos[j{a — 4>j)] + B, where {Aj, cj)j) are the amplitude and phase of the anisotropy, 
a is the right ascension, and B is a constant. Fig. [2] shows the right ascension projection of 
the maps from Fig. [1] for the two energy samples together with the fit result. The error bars 
represent the statistical errors. 

The anisotropy observed at 20 TeV with IC59 is consistent with the previously reported 
observation with IceCube-22 confirming the continuation of the arrival distribution pattern 
observed in the northern hemisphere to the southern sky. At 400 TeV the sky map and the 
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Figure 1: Relative intensity AN/{N) of the IC59 data in equatorial coordinates for the 20 TeV (left) and the 400 
TeV (nght) energy sample. The sky maps are smoothed within 3° to match the angular resolution of the data. 
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Figure 2: One-dimensional projection in right ascension a of the 20 TeV (black) and 400 TeV (red) maps in Fig.[T] 
The data are shown with statistical uncertainties, and the lines correspond to the first and second harmonic fit. 



right ascension projection show substantial differences with respect to the observation at lower 
energy. The large scale structure visible at 20 TeV vanishes and a smaller (of ~ 20° extension) 
deficit appears at (a = 73.1°, S = —25.3°) with a post-trials significance of —6.3a. A significant 
excess cannot be distinguished yet from the current statistics. This is the first observation of an 
anisotropy at 400 TeV in the southern sky^l^. 

3.2 Medium-Small Scale Structures 

The arrival direction distribution of cosmic rays observed in IC59 is dominated by the presence 
of large structures in the form of a strong dipole and quadrupole component. In order to search 
for correlations on smaller angular scales a complementary method has been used to analyse the 
data. The first step in this search is the creation of a "reference map" to which the actual data 
map is compared. The reference map essentially shows what the sky map would look like if the 
cosmic ray fiux was isotropic. The reference map accounts for detector effects that might cause 
a spurious anisotropy, such as non-uniform exposure to different parts of the sky or gaps in the 
uptime, and it is created using a background estimation method based on the time scrambling 
technique described in Alexandreas et al.^^] 

The complete analysis, reported in R. Abbasi et a can be summarised as follows: 

• The muon events are binned in equatorial coordinates to create a 2D map of event counts 
Ni{a, 6). (Note that i is the pixel index in the binned map.) The reference counts {Ni(a, 6)) 
are estimated from the data, and deviations from isotropy are calculated in each bin. The 
deviations can be expressed in terms of a relative intensity 6Ii = (Ni — {Ni))/{Ni), or in 
terms of statistical significance using the method of Li and Ma^^. 

• The angular power spectrum of the relative intensity map is calculated to determine corre- 
lations between pixels at several angular scales. The angular power spectrum of the IC59 



data exhibits significant power not only at the largest angular scales, but also down to 
scales of about 10°. 



• A fit of the dipole and quadrupole terms to the relative intensity map is performed and 
subtracted from the intensity map in order to obtain a map of residual counts, which is 
then analyzed for small-scale structure. 

• To increase the sensitivity to the small-scale structures in the data, we apply a smoothing 
procedure which takes the reference level and residual data counts in each bin and adds 
the counts from pixels within some angular radius of the bin. 
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Figure 3: Left: Significances of the IC59 residual map plotted with 12° smoothing. Right: Significances of the 
IC59 residual map plotted with 20° smoothing. The regions with a pre-trial significance larger than ±5a are 

shown in colours. 

Fig. [3] shows the smoothed residual maps of significance in equatorial coordinates. No energy 
selection has been done for this analysis; the maps show the sky at a median energy of 20 TeV 
(Sec. [2]). We identify eight regions, shown in colours, with a pre-trial significance larger than 
±5(7. Since the optimal scales vary from region to region and no single smoothing scale shows 
all regions, we show the maps with two smoothing scales, 12° (left) and 20° (right). The most 
significant excess on the sky is Region 1 at (a = 122.4°, 5 = —47.4°), with a peak significance 
of 7.0a at a smoothing radius of 22° which is reduced to 5.3fT after accounting for trials. 

4 Conclusions 

Using 3.2 X 10^'^ events recorded with the partially deployed IceCube detector between May 2009 
and May 2010, we have found that the arrival direction distribution of cosmic rays at several 
TeV exhibits significant anisotropy on several angular scales. The data are dominated by dipole 
and quadrupole moments, but there is also significant structure on angular scales down to about 
10°. At 400 TeV the large structure seems to vanish and the most significant region is a relative 
deficit at a right ascension where the broad excess dominated at 20 TeV. 
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